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Abstract 

We study a dilute gas with two species of Fermionic atoms of unequal con- 
centrations, interacting via a short-range interaction with one deeply bound 
state. We study the properties of this system under the mean-field approxi- 
mation. We obtain the effective interaction among the fermions and bosons, 
and discuss the collective modes of the system. 
PACS numbers: 03.75.Fi, 67.60.-g 
Recent experimental advances in trapping and cooling alkali atoms have already given 
us condensates of Bosons [1] and degenerate gases of Fermions [2,3]. Much theoretical and 
experimental efforts have now turned to systems with trapped molecules, with interesting 
observations and predictions on the equilibrium and dynamical properties of this system. 
[4-9] These molecules can have their constituent atoms be either Bosons or Fermions. [10] 
The latter case, on which the rest of this paper will concentrate, is also of particularly interest 
when it is considered together with superfluidity arising from Cooper pairing of Fermions. 
Indeed superfluid 3 He and most known superconductors and have their superfluidity derived 
from Cooper pairs. It is also widely believed that Fermionic pairing in diatomic molecules 
and in Cooper pairs are in fact two limits of a continuum. [1 1] Much work has already been 
devoted to the cross-over between these two limits. [12,13] (and references therein) There is 
also a lot of interest in producing pairs (within any part of the above mentioned continuum) 
in trapped atomic Fermions. [10] 

We here would like to further our understanding of systems consisting of "molecules" . We 
shall limit ourselves to "molecules" made up of two distinguishable Fermions, interacting 
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via a short range (<< a below) interaction, in the s-wave state (Therefore by 'molecule' 
we simply mean the Boson made up by a pair of Fermions). These two partners of the 
molecule can in general be atoms in different hyperfine spin projections or internal states, 
or different atomic species. For definiteness we shall restrict ourselves to the first case and 
for simplicity refer them hereafter as spin "up" and "down", with concentrations n<\ and 
n j. We shall study the properties of this system when the two constituents of the molecule 
are of different concentrations, say when n-f > nj_. Furthermore, we restrict ourselves to the 
case where the size of the bound state a is small, such that n^a 3 and n^a 3 are both much 
smaller than 1. Then at temperatures much below the dissociation temperature Tdi SSO all 
down spins are paired into molecules but there are left over (or excess) spin-up Fermions. 
The Bosonic molecules and Fermionic atoms basically do not overlap and are well-defined 
entities. Thus one has a Bose-Fermi mixture where the Fermions are identical with one of 
the constituent particles of the Bosons, with the Boson and Fermion concentrations rib = n± 
and v = n-f — respectively. 

We shall obtain the properties of our Bose-Fermi mixture by applying many-body tech- 
niques on the (more basic) constituent spin- up and down Fermions (generalizing ref [11-13] 
- we shall also be following the notations of the last reference quite closely, except that 
we shall use the more common sign conventions for the single particle propagators and self- 
energies). This system is attractive theoretically in that all interactions can be characterized 
by one interaction parameter, namely the s-wave scattering length a (> 0) between the spin 
up and down Fermions, or equivalently the binding energy % of the molecule (e& = h 2 /ma 2 , 
here m is the mass of either of the Fermions). There are, however, several energy scales 
and hence temperature regimes. As mentioned we shall limit ourselves to where the tem- 
perature T of the system is much lower than the dissociation temperature Tdi SSO (~ up 
to some entropic factors). The relevant scales are thus T, the Bose-Einstein condensation 
temperature T B ec for the Bosons and the Fermi temperature T F for the (left-over) Fermions 
( T BEC ~ nl^ 3 h 2 /m and T F ~ u 2 ^h 2 /m. Note they are much less than T disso ). We shall 
study the some properties of this system at T = and compare them with those of a general 
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dilute Bose-Fermi mixture. In particular, we shall be able to find the effective interaction, 
or equivalently the scattering length, between the Bosons and Fermions in terms of a. We 
further show that the Bogoliubov mode of the system is the same as that of a Bose-Fermi 
mixture. 

In the case where the pairing between the spin-up and down particles are in the Cooper 
limit (weak attractive interaction so that a < and small), the system with raj > has 
been studied before in the context of superconductivity under an applied Zeeman field h 
[14] . We shall compare our results with those of this better known system. 

We now consider the system at T = under the mean-field approximation The necessary 
mean field equations can be obtained by generalizing those of ref [11] (which deals with 
n T = n l)- We introduce three fields, the chemical potentials /if, ji^ for the two spin species 
and a pairing field A. These fields obey the self-consistent equations 
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[l-f(E k -h)-f(E k + h)) 



(1) 



n a = I n a (k) a = t or J, (2) 
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with 



n } {k) = u\ f(E k - h) + v 2 k (1 - f(E k + h)) (3) 

ni (k) = u\ f(E k + h) +v 2 k (l- f(E k - h)) (4) 

where / is the Fermi function. We have defined for convenience the chemical potential 
average and difference /i = (/i^ + AH)/2, h = — //j.)/2 (> 0), the energies e k = % 2 h 2 /2m, 
E k = [{e k - fi) 2 + |A| 2 ] 1/2 , and the "coherence factors" u\ = |(1 + ^), v 2 = \{l - ^), 
u k v k = are exactly as in BCS theory for superconductivity. Eq (1) is the self-consistent 
equation for A whereas eqs (2) are the constraints for the Fermion concentrations, with 
n a {k) the occupation number at wavevector k. Notice that if we view the system as a Bose- 
Fermi mixture, the chemical potential for the Bosons is /if, = 2/z and that for the (left-over) 



Fermions is /if. In the dilute limit < Xf = r^a 3 << 1 and < x^ = n^a 3 << 1 which we 
are studying (in contrast to the Cooper limit where a < ), it is convenient to study the 
solution to eqs (l)-(2) as an expansion in xj and x^. At T = 0, f(Ek + h) — and /(-£*; — /i) 
is non-vanishing only for k < k c where E^ c = h. From eq (3) and (4) we get immediately 

v = k*/6n 2 . (5) 

k c plays the role of the Fermi wavevector for the excess spin-up Fermions (which form a 
'Fermi sea', see also below) 

We shall see (as in the corresponding case of ref [11]) that ji < and |A| << Using 
this inequality, eq(l) and eq(4), we obtain two equations determinating A and /i. In the 
dilute limit we get 

_ h 2 + 2ir(n b + 2v)h 2 a 
2ma 2 m 

|A| 2 = 8im b e 2 b a 3 [1 - ...] (7) 

(and hence |A| 2 //i 2 ~ n^a 3 << 1 as promised). Eq (6) can be interpreted easily by recalling 
\i h = 2/1. The first term arises from the binding energy of the molecule. The two parts of 
the second term represent the mean-field corrections to the chemical potential due to the 
interaction of the molecules among themselves and with the Fermions respectively. Equating 
these two expressions with (1/2 times) n^bb an d vg b f, we obtain the interaction constants 
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or equivalently the scattering lengths = 2a and a b f = 8a/3. [15] 

We now discuss some properties associated with the Fermions. The occupation numbers 
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of the constituent Fermions n a (k) are sketched in Fig. 1 for T = 0. n^{k) = 1 whereas 
ni(k) = for k < k c . n^{k) and n±(k) are equal but small for k > k c . These results can be 
understood by considering the wavefunction of the sytem [16]. 



It is also interesting to study the spectral function A^(k,uj) (= — ^ Im Gf u (k,u) given 
below) corresponding to the density of states for addition/removal of a spin- up Fermion 
at momentum k and energy uj. A^{k,u) has the same formal expression as in the case of 
Cooper pairing (c.f. eq (12) below): A(k,u) = u 2 k 5{uj — (E k — h)) + v k 5{uj + (E k + h)). 
There are thus two delta-function peaks, one at positive and one at negative energies. Since 
|A| << /i, v\ « 1 and thus the weight at the negative frequency is very small. To a very 
good approximation, ~ 5{uj — (e k — /i|)), that of a free Fermi gas. We shall compare 
these results with those of case (B) below. 

We next study the dynamic properties of our system (at T = 0) by evaluating the two 
particle vertex function (matrix) T(Q, Q) for repeated scattering between a spin- up and 
spin-down particle (hole) at small wavevector Q and frequency Q. We shall see that, as in 
the case of — n±, [12,13] T is related to the propagator for the Bosons. T, in Matsubara 
frequencies (we shall drop the superscript M for simplicity. The desired T in real frequency 
can be obtained via iQ u — > Q with Q u > as usual) obeys the Bethe-Salpeter equation 
(generalizing that of [13]) 

KlAQ, ^v) = -r-T2-S*i,a2 + M alt0a (Q, n u ) (10) 
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and G a a i ta2 are the matrix Green's function in particle (a = 1) and hole (a = 2) space for 
spin-cr defined in a similar manner as in Gorkov's theory of superconductivity [17], and fl u 
and e n are Bosonic and Fermionic Matsubara frequencies respectively. In principle there is 
a contribution to the self-energy for G proportional to T and therefore G and T have to be 
found self-consistently (see [13] ). However, at T = and under our diluteness condition we 
can simply use the mean-field solution for G, i.e., 

G m(k,e n ) = , % + . ft, (12) 

ie n - (E k - h) ie n + (E k + h) 



and 



G n2 (k,e n ) = u k v k 
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(13) 



_ie n -(E k -h) ie n + (E k + h)_ 

where u k , v k are the coherence factors introduced before. Gj/s are given by similar expres- 
sions except h — > —h and A — > —A. r _1 (<5 = 0, VL U = 0) can be calculated easily and is of 
order A 2 due to eq (1). However, evaluation of T _1 at finite Q and ft needs special care. In 
Mn there arises the contribution 
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The integrand is singular as Q and ft approaches 0. To evalute M{\, we first note that under 
our diluteness approximation, we can replace E k by e k — fi and thus E k — h = E k — E kc m 
e fe — £fcc etc. The integrand is therefore finite only if 6q_£ c — e fcc — ft and e k — e kr are of 
opposite signs, and thus for Q « k c , ft << k 2 c /m, k is restricted to be near the "Fermi 
surface", i.e., \k\ ~ k c . We can then use « u\ c v kc « , and replace the integral 

by integration over fc and ^ k = e k — e kc . We finally obtain M{\(<5, ft) ~ f^('~) 2 ! an d hence 
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As in the case of equal number of spin-up and down particles, [13] T is thus found to be 
proportional to the connected part of the Bosonic propagator matrix G&. G^ 1 , under the 
random phase approximation, reads (c.f. [18]) 



-G bl \(Q, ft) — —ft + — + g bb n b - g bf n b Xf{Q, ft) 



(16) 



-G b n(Q, ft) = +g bb n b - g bf n b x f (Q, ft) 
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where Xf is the density response of the Fermions and we have taken the gauge where the 
phase of the condensate is zero. To further check our claim that r oc G b above, we note that 
m b = 2m, g bb n b = g 2 bf n b = 87r|A| 2 a 3 (using eqs (7), (8) and (9)), Xf{Q,ty = (0). 
In this last expression, we have used the value of Xf(Q,ty f° r &/Q ~^ 0- That this is the 
correct limit can be seen as follows. The bare sound velocity for the Bosons is given by c° = 
^b9bbl m b) x ^ 2 whereas the bare Fermi velocity is v® = k c /m. Thus (c^/v®) 2 ~ (n b a/k 2 ) « 1 
under our diluteness condition. For the mode under consideration, (Q/vfQ) ~ (c b /v j) « 1. 

In conclusion, we have studied, via many-body technique, some properties of an atom- 
molecule mixture where the Fermionic atoms are identical with one of the constituent atoms 
of the diatomic molecules. If T << T BEC) our system behaves like a general Bose- Fermi 
mixture. In this case we have also obtained the scattering lengths a bb and a b f for our mixture 
in terms of a. 

I would like to thank Darwin Chang for pointing out and making available to me ref 
[13]. This research was supported by the National Science Council of Taiwan under grant 
number 89-21 12-M-001-061. 
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FIG. 1. Schemetic plot of n a (k) versus k. Full line: nf(fc), dotted line: n^k). The magnitude 
of n a (k) for k > k c has been exaggerated for clarity. Inset: Sketch of n(k) for the case of nt = n\. 
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